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Adiabatic Continuation between 
Resonating- Valence-Bond Electron and Spin Liquids 

Fusayoshi J. Ohkawa*^ 
Division of Physics, Hokkaido University, Sapporo 060-0810, Japan 

The Hubbard model in the strong-couphng regime is studied by the Kondo-lattice theory. If no 
- - - symmetry is broken at a sufficiently low temperature in sufficiently low dimensions, an electron liquid 

is stabilized by the Fock-type exchange effect of the superexchange interaction. Since the stabilization 
• mechanism is none other than the resonating-valence-bond (RVB) mechanism, the electron liquid is 

■ none other than an RVB electron liquid. The RVB electron liquid in the Hubbard model is adiabatically 
connected to the RVB spin liquid in the Heisenberg model. 
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§1. Introduction 

m 

The concept of adiabatic or analytic continuation is the counter concept of symmetry break- 
^ ■ ing: If no symmetry is broken or restored when a phase transforms into another phase as a function 

of an adiabatic parameter or the temperature T, the two phases are adiabatically or analytically 
\ continued or connected to each other. The concept is useful to understand, for example, the Kondo 

effect. 

The s-d model is derived from the Anderson model in the s-d limit. Since the s-d exchange 
^ _ interaction is antiferromagnetic, the ground state of the s-d model is a singlet because of the 

■ quenching of the localized spin by conduction electrons.'^ '^''^ The Kondo temperature Tk or 
C . /^bTk is the energy scale of local quantum spin fluctuations. The phase transformation between a 

local-moment state at T S> Tk and a local spin liquid at T ^ Tk is a crossover, and no symmetry 
is broken in the lower temperature phase. Thus, the local-moment state and the local spin liquid 
are analytically connected to each other.EJ Furthermore, the local spin liquid can be described as 
the local normal Fermi liquid.'^ Since local charges carried by the localized spin cannot fluctuate 
in the least, the existence of charge fluctuations is not a necessary condition for the normal Fermi 
^sgj . liquid to be stabilized. 

OO \ Since electron correlations are local in the Anderson model, no cooperative effect is possible 

^ ■ and no symmetry can be broken in it; therefore, the conventional perturbative analysis in terms 

of the onsite repulsion U is valid, no anomalous term appears in the perturbation, and the ground 
^ \ state is the local normal Fermi liquid for any finite ?7 > 0, even in the s-d Thus, the 

local Fermi liquid for any finite C/ > is adiabatically connected to that for [/ = 0. The Fermi- 
liquid relation for the Anderson model in the s-d limit is exactly the same as that for the s-d 
^ I model.^'-'*''^''^'''' The facts discussed above mean that the normal Fermi liquid in the Anderson 

^ ■ model for any [/ > is adiabatically connected to the spin Uquid in the s-d model, even if local 

charges can fluctuate in the Fermi liquid in the Anderson model but they cannot, in the least, in 
the spin liquid in the s-d model. On the basis of the adiabatic continuation, it is reasonable that 
physical properties resemble each other between the Anderson and s-d models. 

The relation between the Hubbard and Heisenberg models is similar to that between the An- 
derson and s-d models. If C//|i| » 1 in the half-filled Hubbard model, where U is the onsite 
repulsion and — f is the transfer integral between nearest neighbors, the band splits into the upper 
and lower Hubbard bands, which are empty and full, respectively;^''^' and a narrow midband pre- 
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dieted by Gutzwiller's theorjPJ'EJnD ean appear within the Hubbard gap between them.'^ On 
the other hand, almost all unit cells are singly occupied and the probabilities of empty and double 
occupancies are almost vanishing but still nonzero. The superexchange interaction arises from the 
virtual process allowing empty and double occupancies.^"*' In field theory, on the other hand, any 
mutual interaction arises from the virtual exchange of a boson or a bosonic excitation; e.g., the 
nuclear force arises from that of a Yukawa's meson or a pion. The virtual process from which the 
superexchange interaction arises is none other than the virtual exchange of a pair excitation of an 
electron in the upper Hubbard band and a hole in the lower Hubbard band. The pair excitation is a 
bosonic excitation. The superexchange interaction arises from the virtual exchange of the bosonic 
excitation. The superexchange interaction is antiferromagnetic and the exchange interaction 
constant between nearest-neighbor spins is J = — 4t^/[/. Thus, the Heisenberg model is derived 
from the Hubbard model in the Heisenberg limit. It is anticipated that physical properties resemble 
each other between the half-filled Hubbard model with U /\t\ » 1 and the Heisenberg model with 
J = —At^/U, except for physical properties directly related to charge fluctuations such as the 
electrical conductivity. 

The role of the superexchange interaction is dual: the cause and suppression of magnetic 
instability. Since the superexchange interaction is antiferromagnetic, it is possible for an antifer- 
romagnetic state to be stabilized below the Neel temperature T^. On the other hand, Fazekas and 
Anderson proposed the resonating-valence-bond (RVB) theory for the Heisenberg model on the 
triangular lattice. A spin liquid is stabilized at least at T > OK because of the formation of 
resonating valence bonds between nearest-neighbor spins by the superexchange interaction. The 
spin liquid is the RVB spin liquid and the stabilization mechanism is the RVB mechanism. 

The RVB mechanism is also effective in not only the triangular lattice but also other types 
of lattice in sufficiently low dimensions. The energy scale of resonating valence bonds or local 
quantum spin fluctuations on nearest neighbors is 0(| J|). If Tn ^ \J\/k^, a spin state at a 
temperature T such that Tn < T <^ | J|/fcB is a spin liquid stabilized by the RVB mechanism, i.e., 
the RVB spin liquid. Thus, it is anticipated that an electron state at the T in the Hubbard model that 
corresponds to the Heisenberg model is an electron liquid stabilized by the RVB mechanism, i.e., 
the RVB electron liquid.^"^ The RVB mechanism plays a crucial role in quenching magnetic 
moments, as the Kondo effect. It is interesting to examine whether the RVB electron and spin 
liquids in the Hubbard and Heisenberg models are adiabatically connected to each other, as the 
local electron and spin liquids in the Anderson and s-d models. 

On the other hand, the Kondo-lattice theory has been developed in previous paperi^''^''^''^''23''S) 
in order to study strongly correlated electrons in lattice models where the onsite U is crucial, 
such as the periodic Anderson model and the Hubbard model. For example, the self-energy of 
the Hubbard model is decomposed into the single-site and multisite self-energies; the single-site 
self-energy can be mapped to the local self-energy of the Anderson model; then, the multisite 
self-energy can be perturbatively calculated on the basis of the local self-energy of the mapped 
Anderson model. On the basis of the mapping to the Anderson model, the Kondo temperature 
Tk or /cbTk can be defined as the energy scale of local quantum spin fluctuations in the Hubbard 
model. The k^TK is also the energy scale of the effective Fermi energy; if T ^ Tk then electrons 
behave as itinerant electrons while if T » Tk then electrons behave as local moments. Thus, if 
Tn ^ Tk in sufficiently low dimensions and if Tn < T ^ Tk, an electron liquid is stabilized by 
the Kondo effect. The Kondo effect also plays a crucial role in suppressing magnetic instability 
in lattice models where the onsite U is crucial, as the RVB mechanism. It is interesting to study 
by the Kondo-lattice theory how the Kondo effect and the RVB mechanism cooperate with each 
other to suppress magnetic instability. 

One of the purposes of this paper is to study the RVB electron liquid in the Hubbard model. 
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The other purpose is to show that the RVB electron liquid in the Hubbard model is adiabatically 
connected to the RVB spin liquid in the Heisenberg model. This paper is organized as follows: 
The Hubbard model studied in this paper is defined in f|2] In f|3l the Kondo-lattice theory is 
reformulated in order to treat properly the RVB mechanism. In ^ it is shown that the Kondo 
temperature is enhanced by the RVB mechanism, so that k^TK oc \J\ and fesTk = 0(| J|) for 
the RVB electron liquid. The adiabatic continuation between the RVB electron and spin liquids 
is examined in ^ Discussion is given in ^ Conclusion is given in ^ A sum rule is proved in 
AppendixEl Two applications of the sum rule are given in Appendix|B] An equality is proved in 
AppendixICl 

§2. Hubbard Model 

The Hubbard model on a bi-partite hypercubic lattice in D dimensions is studied in this paper: 

•i(T ija i 

where and d^^ are creation and annihilation operators of an electron at the ith. unit cell, riifj = 
dl^d-^, ed is the band center and = is assumed for simplicity, —t/VD is the transfer integral 
between nearest neighbors, 6(^ij^ = 1 between nearest-neighbor unit cells and 6(^ij) = between 
non-nearest-neighbor unit cells, U is the onsite repulsion, /x is the chemical potential, and 



M = Y,^i.- (2) 



la 



The periodic boundary condition is assumed. If U = 0, the dispersion relation of an electron as a 
function of wave number fc = (fci, /c2, • • • , A:^) is given by 

1 ^ 

E{k) = ed- 2t(pD{k), (foik) = —= ^ zo^{Ka), (3) 



where = 0, as discussed above, and a is the lattice constant. The density of states is given by 

1 

I 



Po{e) = j^S[e + fi-E{k)], (4) 



k 



where L is the number of unit cells. The definition of the transfer integral in Eq. ([T]) is different 
from that in fJTJ it includes the dimensional factor of 1/ as in previous papers.'^''^''23JIl} 
Since (1/L) J2k "Poi^) = the effective bandwidth of £;(fc) or po(e) is 0{\t\) for any D. 
The electron density per unit cell is defined by 

n = {M)/L, (5) 

where (• • • ) stands for the statistical average. If ^ = + U /2, then n = 1, i.e., the Hubbard 
model is half filled; and there exists the particle-hole symmetry in it, i.e., the Hubbard model is 
symmetrical. It is assumed in this paper that ^ is one such that n = 1, n = 1 it 0+, or n ~ 1 for 
it, and that U /\t\ is in the strong-coupling regime defined by 

U/\t\ > 1. (6) 



The thermodynamic limit of L — )• +oo is assumed. No symmetry can be broken at nonzero 
temperature in one and two dimensions,"^ and symmetry can be broken below a nonzero critical 
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temperature Tc in three dimensions and higher. It is assumed that the dimensionality D is suffi- 
ciently low such that Tc does not exists or, if it exists, <^ Tk, where Tk or /cb^k is the energy 
scale discussed in fJT] and the definition of it is given in ^3.1[ the case of Tc ~ Tk or Tc > Tk is 
out of scope in this paper. Then, it is assumed that the temperature T is within the range of 

Tc + 6T<T<Tk, (7) 

where 6T is sufficiently small such that 6T <^ Tk, but still sufficiently lai^ge; e.g., 6T ~ O.ITk, 
or 5T ~ O.OITk, or 5T ~ O.OOITk, depending on D and /U, or n0 No symmetry can be broken 
in the temperature range defined by Eq. Q. 

If no symmetry is broken, the Green function in the wave-number representation is given by 

Ga {iei , k) = - — ■ — — — \—. — — — — , (8) 

lei + ii- E{k) - E„{iei,k) - r{iei) 

where ei = {21 + l)TTk-QT, with / being an integer, is a fermionic energy, and Efj{iei, k) is the 
self-energy. If /i = + [7/2 in two dimensions, i.e., if n = 1 and D = 2, po{e) diverges 
logarithmically as e — 0. In order to suppress the divergence, r{iei) is introduced; it is due to, 
e.g., an explicit electron reservoir such as one considered in the previous papei^ or impurities 
considered in ^4.41 of the present paper. For the moment, it is assumed that r{iei) is almost 
infinitesimal but still nonzero; i.e., \mr{e + iO) = —0+ for e ~ 0, where r{e + iO) is the retarded 
one of r{iei). 

The site-diagonal Green function, the density of states, and the electron density are given by 

Y r+oo 

RAi^i) = T G^{iei,k), p{e) = — Im/?,(e + iO), n = 2 def{e)p{e), (9) 

-^-^ 

respectively, where {e + iO) is the retarded one of i?cr ) and 

f{e) = l/[e'l^^^'^^ + 1]. (10) 



§3. Kondo-lattice theory 



3.1. Mapping to the Anderson model 
3.1.1. Self-energy 

The decomposition of the self-energy into the single-site and multisite ones can also be de- 
scribed in the wave-number representation: Ua-ii^i, ^) = ^aii^i) + ^^a{i£u where Ufj{i£i) 
is the single-site one and AIlfj{iei,k) is the multisite one. The single-site self-energy can be 
mapped to the local self-energy of the Anderson model.E^''^''22 The Anderson model to be 
mapped or the mapped Anderson model is specified by four parameters of T, U, Id — p, and A{£), 
where T is the temperature of the reservoir for the Anderson model, IJ is the onsite repulsion, 
which is denoted by ?7 in 21 is the level of localized electrons, ji is the chemical potential, and 
A{e) is the hybridization energy between localized and conduction electrons. The Green function, 
the density of states, and the density of localized electrons are given by 

-1 

(11) 



*' If Tc does not exists, Tc — Q has to be assumed in Eq. l[7]l. The ST is introduced partly because comphcation 
arising from the logarithmic divergence of po(e) as e — >■ can be avoided and partly because critical phenomena are 
out of scope in this paper, except for an argument in ij6] in which a possible anomaly in the critical region in two 
dimensions is examined. 
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p{e) = — lmG„{e + iQ), n = 2 def{e)p{e), (12) 

71" J-oo 

respectively, where ei = {21 + l)iTkBT, with I being an integer, and ZJa-ii^i) is the local self- 
energy. According to previous papers,?^''-"*' if the four parameters are given by 

f = T, il = U, id - p = ed - p, (13a) 
A{e) = Im [U^{e + iO) + + iO)] , (13b) 

it follows that iei = iei for any /, and 

Kiiei) = Kiiei), (14a) 
G^{iei) = Ra{iei), p{e) = p{e), h = n. (14b) 

Not only the single-site Uo-{i£i) but also the local Ra-{i£i), /o(e), and n are mapped to the local 
Ila{iei), Gjj{iei), p{e), and n, respectively, of the Anderson model. Equation ([T3] ) is the mapping 
conditionO Since A{£) determined by Eq. (I13bl ) depends on the temperature T of the reservoir 
for the Hubbard model, the mapped Anderson model depends on the T, i.e., the Anderson model 
includes the T as a parameter. If the Hubbard model is half filled and symmetrical, the Ander- 
son model is also half filled and symmetrical, so that h = 1, A{e) = A{—e), p{e) = p{—£), 
lxnG^{£ + iQ) =ImG'<^(-e + iO),ReG'<^(e + iO) = -ReG'^(-e + iO), ReG^i+iO) = 0, and so 
on. 

If the multisite self-energy is ignored in Eq. (I13bl ). i.e., if AIlfj{iei, fc) = is assumed, the 
Kondo-lattice theory is reduced to the supreme single-site approximation (S^A),'23'I23J22)J2SI which 
is rigorous for infinite dimensions or, precisely speaking, for 1/D = 0]**i but within the con- 
strained Hilbert subspace where no symmetry is allowed to be broken j***1 Either the dynamical 
mean-field theorjE^I^ or the dynamical coherent potential approximation^^ ' is the S'^A theory. 

The mapping condition of Eq. (I13bl ) can also be described as 

A{^) = ^ - |lmZ'<,(e + iO)|. (15) 

[Rei?,(e + fO)]'+ [7rp(e)]' ' ' 

If U /\t\ = 0, + ^0, k) = 0; and either Rei?o-(£ + jO) or po(^) is continuous and finite at 
e = 0, even fox D = 2 and n = 1 in the presence of -r(e + iO). Then, it follows from Eq. ([T5] ) that 

< lim Aie) < +oo, (16) 

e->0 

for [7 = 0. On the other hand, according to the previous paper,!^ 

^(e) > |lmr(e + iO)|, (17) 

is satisfied for any finite U. A critical Uc can be defined such that Eq. ([T6l ) is satisfied for < ?7 < 
Uc while it is not satisfied for U > Uc- In the temperature range defined by Eq. (|7]l, no symmetry 



*' In the actual mapping procedure, A{£) has to be self-consistently determined with the single-site and multisite 
self-energies according to Eq. l ll3bt in order for Eq. ( I14ab to be satisfied. 

**' It is quite easy to show that, under S'^A, if Eq. il6\ is satisfied then [p(0)] j,^qj^ is constant as a function of U/\t\, 
i.e., [p(0)] y^Qj^ — po(0). As studied in i]4.2l of this paper, if the RVB mechanism is considered beyond S'^A and if 
Eq. l ll6b is satisfied, p{0) decreases as U/\t\ increases for any finite D, even for 1/D = 0"''; in the half-filled case, in 
particular, p{0)\t\ as U/\t\ +00 unless 1/D = 0. 

• ••) jjjg conventional Weiss mean fields for spin density wave (SDW) or magnetism, charge density wave (CDW), 
and isotropic s-wave or BCS superconductivity are of the zeroth order in 1 /D. Since the conventional Weiss mean 
fields are multisite effects, they cannot be considered in any single-site approximation. 
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can be broken for any finite U; thus, it is anticipated that the conventional perturbative analysis in 
terms of U is valid for any finite U, and that Eq. ([T6l ) is therefore satisfied for any finite U. For the 
moment, it is assumed that the critical Uc is infinite, even for D = 1 and n = 10 The relevance 
or irrelevance of this assumption is critically examined in ^4.3l and Appendix|Bl but independently 
from different points of view; in particular, that in the case of D = I and n = 1 is critically 
examined in ^ 

3.1.2. Polarization function in the spin channel 

The decomposition of the irreducible polarization function iTg {iuji , q) in the spin channel 
into the single-site and multisite ones can also be described in the wave-number representation: 
7rs{iuji,q) = ■Ks{iijJi) + AiTs{iijJi,q), where Trs{iuJi) is the single-site one, A'irs{iLOi, q) is the 
multisite one, and loi = 21'irkBT, with / being an integer, is a bosonic energy. The single-site 
T^sii^i) is also mapped to the local TTs{iuJi) of the Anderson model: TTs{iuJi) = tXs{wi). The spin 
susceptibilities of the Anderson and Hubbard models are given by 

2Tis{iu}i) 2'Ks{iu)uq) 

Xs{^^^l) = -^ ,. Xs\^^^hq) = -^ j^-r. r, (18) 

\-\Jt^s\i^i) l-UTrs{iuJi,q) 

respectively, where the conventional factor (l/4)g'^/i| is not included. 

The Kondo temperature Tk can also be defined for the Hubbard model by 

[UO;T)]^^^^ = l/[kBTK{T)]. (19) 

In Eq. ([T9l ). the subscript T — K means that the temperature of the reservoir for the Anderson 
model is zero; and the temperature of the reservoir for the Hubbard model, i.e., the parameter T is 
explicitly shown because Xs(0) and Tk depend on it. The k^TY^{T) or feelk defined by Eq. ( [T9l ) 
is the energy scale discussed in ^and ^ 

The analysis so far is valid for any finite U /\t\. The analysis in the following part is only valid 
for U/\t\ > 1. It is obvious that if U/\t\ > 1 then kBTK/U < 1. Then, it follows that 

Xsii^i) = 0[l/{kBTK)], Xsii^uq) = 0[l/{kBTK)l (20) 
for T <Tk and \uji\ < A;bTk; and it follows that 

msim) = 1 + 0{kBTK/U), UTTs{i0Jl,q) = 1 + 0{kBTK/U), (21) 

and, therefore, U'^A7rs{iuJi,q) = 0(1). If Eq. (|2TI) is used, it immediately follows that 

U[l - msiiooi)] = [2/xs{i^i)] [1 + 0{kBTK/U)] , (22) 

Xs{iuJi,q) = n/^^w'"^'^ ^~ f ^ [1 + 0{kBTK/U)] , (23) 

I- (l/4:}I,{^UJl,q)xs{^u}l)^ ^ 

for T < Tk and \ui\ < A^bTk, where 

Is{ioJi,q) = 2U^ATr{iL0i,q). (24) 

The terms of 0{kBT\<i/U) in Eqs. (l22l) and (1231 ) are ignored in this paper. Equation (1231 ) is con- 
sistent with the physical picture for Kondo lattices that local spin fluctuations interact with each 
other with an intersite exchange interaction; Isii^^i, q) is none other than the intersite exchange 
interaction. 



*' If a complete gap opens in the self-consistent p(e), as proposed by Lieb and Wu^^' for the half-filled case in 
one dimension, Eq. ( I16t is not satisfied; according to the previous paper,!^ a necessary condition for a complete gap to 
open is that A{Q) = or zi(e) includes the delta function 5(e). 
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The Neel temperature Tn can be determined by Eq. (1231 ): 

Tn = max[TN(q)] , (25a) 

where Tj^{q) as a function of q is defined by 

[1 - {1/4)1,(0, Q)UO)]T=TMg) = 0. (25b) 

In this paper, it is assumed that D is sufficiently small such that Tn ^ Tk or no Tn exists. 
3.1.3. Three-point vertex function in the spin channel 

The reducible and irreducible three-point vertex functions in the spin channel are also decom- 
posed into the single-site and multisite ones. The single-site functions of them can also be mapped 
to the local vertex functions of the mapped Anderson model, respectively. If they are denoted by 
As{iei,iei + iujir,iu)ii) and Xs{iei,iei + iujii;iujii), respectively, then 

As{iei,i£i + iujir,iu}i') = \s{iei,iei + iujir,iuji,)/[l - UTTsii^u)]- (26) 

If Eq. (l22l) is used, it follows that 

U As {iei , iei + iui' ; iui' ) = [2/ Xs {i^v )] {m , m + i^^v ; ) • (27) 

3.2. Normal Fermi liquid in the mapped Anderson model 

If Eq. (fT6l ) is satisfied, the Fermi surface exists in the conduction band of the mapped Ander- 
son model, as discussed in Appendix lB.2[ so that definitely Tk > K. If Tk > K, the ground 
state of the Anderson model is the local normal Fermi liquid, as discussed in ^ even if the mul- 
tisite self-energy of the Hubbard model is anomalous. 

An infinitesimal external Zeeman energy is introduced into the mapped Anderson model: 
T-Lz = —h{h^ — h^), where h = 0+, and is the number operator for localized electrons with 
spin a in the mapped Anderson model. If Tk > K, the local self-energy of the mapped Anderson 
model is analytic on the real axis, so that it can be expanded in terms of iei P'^tS 

IJ^iiei) = To + (1 - + a{l -^s)h- ^{ieif / [it A{{))\ 

- i{ei/\ei\)[Mieif + h2{kBff]/[iTA{0)] + ■ ■ ■ , (28) 

where 0i > 1, </>s > 1, 4'2 ^ 0, (f)2i > 0, and (p22 > 0, in general; and 02 = in the symmetrical 
model. If U/[ttA{0)] > 1, then 0i > 1, 02i > 1, and (^22 > 1; and 

Ws = ^s/^1 = 0{1), W21 = 4>2l/4>1 = 0{1), W22 = h2/^l = 0{l). (29) 

The Ws is none other than the Wilson ratio. If /1(e) is constant as a function of e, Wg = 2 in the 
s-d model or the s-d limit of the Anderson model.'^''^''^ It is anticipated that if U /\t\ S> 1 and 
n ~ 1 then Wg — 2 for the mapped Anderson model, whose depends on e, in general. 

Since the self-energy ZJaiisi) is analytic on the real axis, the Fermi-liquid relatior|23''S).ll)JS) 
has to be satisfied in a self-consistent solution, even if the multi-site AUa-i^ + ^0, k) is anomalous. 
Here, T is treated as being independent of T; the parameter T is explicitly shown. If T — K is 
assumed for a given T, the susceptibility of the Anderson model is given by 

[Xs{0;T)]f^,^ = 1/[A:bTk(T)] = [2,^,(T)/5(0; T)] j^, (30) 
where Eq. ( fT9l ) is used. If T is sufficiently low such that T ^ Tk, then it follows that 



X.(0;T) = [1-0(4)]/[A:bTk(T)] =20,(T)/5(O;T)[1-O(4)], (31) 
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where 5f = f/TK{T); therefore, it follows that 

p{0;T) ~ l/[2MT)kBTK{T)] oc l/[<^i(r)A:BrK(r)] . (32) 

On the other hand, if n = 1 and T ^ K, then ReG^(+iO;r) = and ImZ'^(+iO; T) 0; 
thus, it follows that according to Eqs. ([T4l ) and ([T5l) that 

[7Tp{0-T)A{0;T)]^^^^ = l. (33) 

Therefore, if f < Tk{T) and n ~ 1 then 7rp(0; T)A{0; T) ~ 1, so that 

7r^(0; T) ~ 2^s{T)kMT) cc ^i{T)kBTK{T) . (34) 

The bosonic energy for T is denoted by Cji = 2TTlk-QT, with / being an integer. According to the 
Ward relation,'29 it follows that 

As{iei,iei + iujv-iujir,T) = 1 - lini ^ [d/{dh)\ [^{iei;T) - IJ^{iei;T)\ , (35) 



where a);/ = 0. According to Eqs. (|28] ) and 

As {iei , iei + wv ; wv ; T) = 0s (T) , (36) 
where uii' = 0. According to Eqs. (1271 ) and (l36l ). 

U\s{ieui£i + i(bir,iu]ii]T) = 2(j)s{T) /xs{ii^i']T), (37) 

where a)// = 0. In this paper, Eq. (|37] ) is used for |q| < k-BT]^{T) and < feB^Kir). 

When the above relations are used in a self-consistent solution for the Hubbard model, T = T 
has to be assumed. Every single-site property depends on the parameter T. If the parameter T is 
much lower than Tk, the parameter-T dependence is so small that it can be ignored, except in the 
case of D = 2 and n ~ 1, as discussed in ^ 

If Eq. (|28] ) is used, the Green function of the Hubbard model is given by 

Ga{i£hk) = {l/(i3i)g^{ieuk), (38a) 

Ofjiieuk) = ^ , (38b) 

iei + fi* -[E{k) + AU,{ieuk)\/<t^i-lK{iei) 

/i* = [^-^o)]M, (38c) 
lK{iei) = -i{ei/\ei\)[W2i{ieif + W22{kBT?]k/[^M^)\- (38d) 

In EQ.(l38b]). r{i£i), h, and -^2{i£if /[t^A{0)] are ignored. The Green function given by Eq. 
(|38] ) can describe only electrons in the vicinity of the chemical potential but it cannot describe 
electrons in the upper and lower Hubbard bands; Eq. (I38bl) is accurate for |e;| <^ k^TK and 
T < Tk and is approximate for < ^bTk or T < Tk- If < 1/(^i < 1 and < /cbTk <. \t\, 
a narrow band appears in the vicinity of the chemical potential. The narrow band is none other 
than the Gutzwiller band.^'^nJ'CS The spectral weight of the Gutzwiller band is l/cpi. 

3.3. Superexchange Interaction 

The intersite exchange Interaction I{iuji,q) can be decomposed into three terms: 

Is (0, q) = Js (0, q) + Jq (0, q) + A{0, q) , (39) 
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where Js{0,q) is the superexchange interaction, and jQ{0,q;T) is an exchange interaction that 
arises from the virtual exchange of a pair excitation of an electron and a hole within the Gutzwiller 
band,® ''''''' Eland A{iuji,q) is the sum of all the remaining terms of Is{iuji, q), such as the mode- 
mode coupling term, terms due to quantum and thermal critical fluctuations, and so on. Only the 
superexchange interaction Js{iuJi, q) is considered in this subsection. 

According to Hubbard's theory,'^''2J the band splits into the upper and lower Hubbard bands. 
Since Hubbard's theory is under the single-site approximation, this result on the band splitting 
can be approximately used for high-energy local properties of the mapped Anderson model. The 
Green function and the self-energy of the mapped Anderson model are approximately given by 

{iei) = -— 1 (40a) 

'"^"-'^^ +. ^"-^ (40b) 



i£i + fi-€d iei+ n - ed-U' 

for \ei\ ^> feB^K, where {fi^) is the number of localized electrons with spin a in the presence of 
the infinitesimal Zeeman energy h; in Eq. (I40bl ). A's in the denominators are ignored because they 
are not crucial. If the rigorous (n^) and (n^) are used in Eq. (l40l ). i.e., if 

XsiO) = lhji{d/dh){{h^) - (n^)), (41) 
is used, it follows according to Eqs. (|27] ). (|35] ). and (l40l ) that 



U\s{iei,iei +iuir,iuv) = - ~ -. — ■ h -. — ■ , (42) 

Gl{i£i) \i£i + fJ'-ed i£i + fj, - ed - U J 

for cj// = and le^l » A;bTk. 

The superexchange interaction arises from the virtual exchange of an electron in the upper 
Hubbard band and a hole in the lower Hubbard band, as discussed in ^ and is a second-order 
effect in —tj^TD. According to Eq. (I24l ). therefore, the superexchange interaction is given by 



Jsii^uq) = \Y. e"'-^'^'-'^'^Jvim), (43a) 



I' 



where Ra-{iei) = Ga{isi). If Eq. (I42l ) is used, the static part of Jij{iiOi) is given by 

J^j{0) = 5^ij^J/D, J = -Af/U. (44) 

This agrees with the one derived by the conventional theory.^ 

Since Jij{u} + iO) — as a; — +oo, and Jij{oj + iO) is analytical in the upper-half complex 
plane, Jij(iuji) can be described, in general, as 



-/ dxXj{x)l-^--. , / dx-^ = -. 

DJq \iuji + x lUJi-xJ Jq X 2 



*' According to a previous paper,^ Jq{0, q) has a novel property such that its strength is proportional to fceTk, 
i.e., the bandwidth of the Gutzwiller band. 
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Since Xj{x) has apeak around x = U,it is assumed in this paper that Xj{x) = {l/2)U5{x — U). 
Then, it follows that 

Js{iuji,q) = ^^D{q)^ l-^—jj - -^—jj] ■ (46) 
In the static limit of |a;/|/[/ — )■ 0, Eq. (I46l ) is reduced to 

Jsiiuji,q) = —=ipD{q) = — ^cos{qi,a). (47) 

i/=i 

The Js{iiOi, q) is of higher order in 1/D for almost all the q's, although it is of the zeroth order in 
for particular q's such as q = (0, 0, • • • , 0), q = (7r/a)(ibl, ±1, • • • , ±1), and so on. 

3.4. Perturbative scheme of the Kondo-lattice theory 

Electrons in the vicinity of interact with each other by the mutual interaction mediated by 
spin fluctuations. The mutual interaction of the first order in the spin-fluctuation mode is given by 

U"^ \s {ieii , ieii + ioji ; iuji ) \s {iein , iei» - iui \-iuJi) [xs {i^i , q) - (^^z )] , (48) 

where ieii and iein are the energies of incoming electrons, and ieii + iuJi and iein — iui are the 
energies of outgoing electrons. Since the single-site part is considered in the mapping to the 
Anderson model, it is subtracted in Eq. (|48] ) in order to avoid double counting. It follows that 

Xs{iu}u q) - Xs{i^i) = {l/^)x^{iuji)ll{iuju q), (49a) 

1 - {l/i)IsiiuJi,q)Xs{iuJi) 

If Eqs. (I37l) and ( |49l ) are used for \uji\ < /cbTk, Eq. (|48] ) is simply described as 4)ll*{iuJi,q). 

If once the mapped Anderson model is solved and the single-site IJ^{isi) and Xs{i^i) are 
given, the multisite AUa-{isi, k) and A7:{iuji, q) can be perturbatively calculated in terms of the 
intersite Is{iuJi, q), although the mapped Anderson model has to be self-consistently solved with 
the multisite terms. Since the single-site terms are given by the local terms of the mapped Ander- 
son model, only multisite terms have to be considered in this perturbative scheme in order to avoid 
double counting. In this perturbative scheme, thus, the intersite Is{i(^i, q) has to be treated as the 
bare intersite exchange interaction, and the single-site (ps has to be treated as the bare vertex func- 
tion in the spin channel; and /* {iuji , q) is the renormalized intersite exchange interaction, which 
is enhanced or screened by intersite spin fluctuations depending on q. 

It is straightforward to show that Is{ioJu q) is of higher order inl/D for almost all the q's 
except for particular q's, as Js{ioJi,q); e.g., the /^(O, Q), where Q is the ordering wave number 
determined by Eq. (1251 ). is of the zeroth order in 1/D and it corresponds to the conventional Weiss 
mean field. The Kondo-lattice theory, which is a perturbative theory in terms of Is{iuJi,q), is none 
other than 1/D expansion theory. 

§4. Resonating- Valence-Bond (RVB) Electron Liquid 

4.1. Fock-type exchange effect of the superexchange interaction 

There are two types of self-energy of the first order in Js{iuJi, q). One is the Hartree-type 
self-energy: 

U^J'\iei) = k^TY,e''^'''*k\js{G,m~o{iev). (50) 
V 
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The other is the Fock-type self-energy: 

AUj^^^\tei,k) = ^Y. 4>l\js{iei-ie,,k-p){a'^'^'.cT"''^)G,,{ie,,p), (51) 

Vpa' 

where a = {(7x,CFy,cjz) is the Pauli matrix. The Hartree-type self-energy is included in the con- 
ventional Hartree term, which is given by 

S^^\iei) = k^TY^e'^^'^^^UR^^iiev), (52) 
I' 

as a part of it. Since the conventional Hartree term is one of the terms for the single-site self- 
energy, which is considered in the mapping to the Anderson model, the Hartree-type self-energy 
has not to be considered in order to avoid double counting. What is considered by the Fock-type 
self-energy is the Fock-type exchange effect of the superexchange interaction, which is none other 
than the RVB mechanism/^' thus, the Fock-type self-energy is called the RVB self-energy. 
If Eqs. (l38l) and (|46]l, and the equahty of 

1 ^ 

2if£){k — p) = —r= [cos(A;jya) cos{pua) + sm{kua) sm{p,^a)] , (53) 

are used, the RVB self-energy is calculated to be 

ASf^^'Xiei^k) = {'i/A)^iW^,{J/D)ED{iei)vD{k), (54) 
where Wg = 4>s/4>i is the Wilson ratio, and 

^uiiei) = j^^Y^oipH ^ n{U)g„{iei + U,p)-n{-U)g„{i£i-U,p) 

--/ def{e)-{-. —FT--. -)lmg^{e + iO,p)\, (55) 

IT J-oo 2\i£i-e + U lei-e-UJ J 

where ga{e + iO, p) is given by Eq. (I38bl ). /(e) is defined by Eq. (ITOl ). and 

n(e) = l/[e^/(^BT) _ ^5^) 

It is easy to show that 



lim SD{iei) = Q, lim AS^]^^'^\ieuk) = ^. (57) 

\ei\/U^+oo |e,|/C/->+oo 

In the large limit ofU /\£i\, Sn^iei) is simply given by 

Sd = j^Y^^Dip) j de/(e)(^--Jlmg,(e + iO,p). (58) 

If U /\t\ » 1 and \£i\/U <^ 1, Eq. (1581 ) can be used for S£,{i£i) with a sufficient accuracy. Then, 

AU^J'''^\iei,k) = M^/A)W^{J/D){t/\t\)\SD\^D{k). (59) 

Even if Eq. (|47] ) is used instead of Eq. (|46] ). Eq. (l59l ) can also be derived. The RVB self-energy is 
of higher order in 1/D, as anticipated. 
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Since AU^J^^^\iei, k) given by Eq. (l59l) does not depend on isi, ga{i£i, k) is simply given 

by 

ga{iei,k) = - ■ — , ,. (60a) 

ak) = -2t*ipD{k), 2t* = 2it/4^i) + it/\t\)cji\J\/D), cj = (3/4) (60b) 
where fi* is defined by Eq. (I38cl ) and | J| = At^ /U. The density of states is given by 

p{e) = p{e) = I™ [(lM)5a(e + iO, k) 



irL 

k 



(61) 



Either Eq. (l60l l or (|6TI ) can describe only the Gutzwiller band but cannot describe the upper and 
lower Hubbard bands. 

Since cj or p* , t* , and 0i depend on each other, they have to be self-consistently 

calculated with each other as a funciton of T and p, or T and n(/u). If no symmetry is broken even 
at T = K, the Fermi-surface sum ruld^D-ESI ^^j^ to determine /i* for T = K: 

^(Z^) = Z E / def{e)5[e + p* - m] . (62) 

If r = K is assumed in Eq. (I62l ). Eq. (I62l ) is none other than the Fermi-surface sum rule; e.g., 
certainly = for n = 1, which is required to be satisfied because of the particle-hole symmetry. 
If T = K and //* = are assumed, it is easy to calculate So defined by Eq. (I58l) : e.g., = 
l/n = 0.31831-- - , I ="2! = 2^/2/7r2 = 0.28658 and |='oo| = l/(2i/¥) = 0.283095-- - ; 
i.e., ^1/3 for any D. If = 1/3 and Wg = 2 are assumed, then cj = 1. Since oc 
feB^k and \t*\ = ©(^bTk), as shown in Eq. (|68] ) in ^4.2[ Eq. (I62l ) can be used, approximately but 
with a sufficient accuracy, for T in the range defined by Eq. (O. If n ~ 1 and T is in the range 
defined by Eq. dT]), then <C so that \Sd\ ~ 1/3 and cj ~ 1 or cj = 0(1). 

4.2. Nature of the RVB electron liquid 

The single-site self-energy has to be self-consistently calculated or determined with the RVB 
self-energy, although the self-cosnsistency is not completed in ^4.11 In this subsection, the na- 
ture of the self-consistent solution is studied under the assumption that the self-consistentcy is 
completed. 

If the self-consistent cj)\ is determined for a U such that U/\t\ » 1, then, according to Eq. 



(I60bl) . 

= |t|[(l/,^i) + 2cj|t|/(Z)C/)], (63) 
where cj = 0(1). If U /\t\ » 1 and if n = 1 orn ~ 1, then l/<j)i <^ 1. Thus, it follows that 

{l/4>i) + 2cj\t\/{DU) « 1, 1^1 < (64) 

Since |t*| <C the Gutzwiller band is certainly formed. Its bandwidth is 0(|t*|). 
If U/\t\ » 1 and 1/(^1 is one such that it satisfies 

2cj\t\/iDU) < 1/01 < 1, (65) 

the RVB self-energy is not the main term of Eq. (|63] ). The bandwidth of the Gutzwiller band, 
which is 0(|t*|), decreases as l/(f)i decreases. Ifl/cpi is nonzero in the limit of U/\t\ — )• +00, as 
shown in Eq. (1721 ) in ^4.31 the spectral weight of the Guztwiller band, l/0i , is also nonzero even 
in the limit of U /\t\ — ;> +00. 
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On the other hand, if C//|f| » 1 and \/4>i is so small that it satisfies 

IM < 2cj\t\/{DU) « 1, (66) 

the RVB self-energy is the main term of Eq. (|63] ). Since the electron liquid in the Gutzwiller 
band is mainly stabilized by the RVB mechanism, it is none other than the RVB electron liquid. 
The bandwidth of the Gutzwiller band or the energy scale of the RVB electron liquid is almost 
independent of and is 0\t^ / {DU)\ . However, the spectral weight of the Gutzwiller band, 
1/01, is quite small or almost vanishing. Thus, the RVB electron liquid is almost a spin liquid or 
a quasi spin liquid. 

In the temperature range defined by Eq. ([7]), the Fermi-liquid relations discussed in ^3.2l are 
satisfied, approximately but with a sufficient accuracy for the analysis of this paper. If Eq. (|65] ) or 
(l66l ) is satisfied, it follows according to Eqs. (I60l ) and (|6TI ) that 

p(0) oc l/(,^i|r|), p(O) = O[l/(0i|r|)], (67) 

and it follows according to Eqs. (l32l) and (l67l) that 

A^BTKOciri, A;BrK = o(|r|). (68) 

In particular, if Eq. (I66l ) is satisfied, it follows according to Eqs. (I44l ). (|63] ). and (|67] ) that 

/5(0) oc (Z)[/)/(,^it2) ocD/(0i|J|), p{Q) = 0[{DU)/{^if)\ =0[D/{^i\J\)], (69) 

and it follows according to Eqs. (041), (l63]l, and ([681) that 

^bTk oc t'/iDU) oc \J\/D, kBTK = 0[ty{DU)] = 0{\J\/D). (70) 

The k^Tx is enhanced by the RVB mechanism. The energy scale fceTk of the Kondo-lattice 
theory is essentially the same one as the energy scale | J|/L> of the RVB theory. 

4.3. On the critical Uc and the self-consistent 

It is a crucial issue whether the self-consistent (j)i can be defined in a self-consistent solution. 
This issue is related with another crucial issue whether or not the critical Uc, which is defined in 
a way such that Eq. (fT6] ) is not satisfied for U > Uc, exists and, if it exists, whether it is finite or 
infinite. 

In the Anderson model, in general, if Eq. (fT6l ) is satisfied, then Tk > K; if Eq. ([T6l ) is not 
satisfied, then Tk = K. In general, if Tk > K, the local self-energy can be expanded as in Eq. 
(|28] ) and the expansion coefficient (pi can be defined; if Tk = K, the 0i cannot be defined. The 
contrapositions of these propositions are as follows: If the cpi can be defined, then Tk > K; if 
the (pi cannot be defined, then Tk = K. In general, if Tk = K, the low-temperature entropy is 
so anomalous that the residual entropy is nonzero or the entropy decreases more slowly than the 
T-linear dependence as T — > K. All of these are also true for the mapped Anderson model. 

If the critical Uc exists and it is finite, the self-consistent (pi cannot be defined for U > Uc, 
so that Tk = K for ?7 > Uc- If Tk = K, the low-temperature entropy has to be anomalous, 
as discussed above. Since it is unlikely that a phase with the anomalous low-temperature entropy 
is stable for a finite U such that Uc < U < +oo, it is unlikely that the critical Uc is finite. This 
argument implies that if the critical Uc exists for n = 1 or n 7^ 1 then the critical Uc has to be 
infinite. 

If the RVB self-energy is explicitly considered as in this paper, in particular, it is impossible 
that the critical Uc appears in the range of |t| <C t/c < +00, as discussed below. If the critical Uc 
is finite and Uc ^ \t\,it follows that 

[kBTK]^^^ _^+^\t*\ = i\t\/^i) + 2cjty{DUc), [A;BTK]^^f, ,0+ =0. (71) 



14 



Unless DUc is infinite, even if a.s U ^ Uc — 0"'", k^Tx is discontinuous at the critical 

Uc', tfius, tfie transition at tfie Uc is a first-order or discontinuous transition. The nonzero k^Tx 
means that the energy gain due to the Kondo effect assisted by the RVB mechanism is nonzero, 
while the zero /cbTk means that the energy gain is zero. Therforer, it is anticipated that the free 
energy of the phase at U = Uc — whose /cb^k is nonzero, is lower than that of the phase at 
U = C/c + 0^, whose k^Tx is zero. If this anticipation is true, it contradicts the possibility of 
the discontinuous transition, so that it contradicts the assumption that the finite Uc exists. If the 
critical Uc appears in the range of f/c 3> ItL the critical Uc has to be infinite unless l/D = 0. 
According to Gutzwiller's theory,E^'[n)'[n) which is for the canonical ensemble, 

lim 1 /4>i = Cc I (N/L) - 1 1 , (72a) 

U/\t\^+oo 

where Cc > and Cc = and is the number of electrons and is an integer. Since 

Gutzwiller's theory is under the single-site approximation, the RVB mechanism is not consid- 
ered in it. Either if the RVB mechanism is considered or if it is not considered, it is anticipated 
that at least 

1/(^1 > Cg|n- l|, (72b) 

where Cg > and Cg ~ has to be satisfied for finite or infinite U/\t\ in the grand canonical 
ensemble. Thus, it is certain that the self-consistent l/4>i can be defined at least for n ^ 1, even 
for n = 1 lb 0+, and for finite or infinite U/\t\. If n / 1, the critical Uc does not exist. 

According to the analysis in Appendix lB.2[ if the RVB mechanism is considered and if n = 1 
orn = 1 lb 0+, then 

lim l/4>i = 0[t^/{DU^)]. (73) 

U/\t\^+oo 

If n = 1 orn = 1 lb 0"'", the critical Uc is infinite. 

On there other hand, if U/\t\ » 1 and n = 1, the probability of empty or double occupancy 
is O [t^/ (DC/^)] . If U /\t\ 3> 1 and n is very close to unity, it is anticipated that the probability is 
also O \p' I (Df/^)] . If U l\t\ ^> 1 and n is sufficiently different from unity but is still sufficiently 
close to unity, the probability is proportional to |n — 1|. Thus, it is anticipated that the probability 
is as large as max[cg|n — 1|, C0t^/(DC/^)] , where c,^ > and c<^ = 0(1). It is a reasonable 
conjecture on the basis of Eqs. d72l ) and (173] ) that the spectral weight of the Gutzwiller band, l/0i , 
is proportional to the probability of empty or double occupancy: 

1/(^1 oc max[cg|n- l|,c^tV(^f^^)]- (74) 

If this conjecture is true, the electron liquid within the Gutzwiller band is the RVB electron liquid 
provided that 

|t|/?7«l, \n-\\<\t\l[DU), (75) 

are satisfied. The RVB electron liquid can be stabilized only in the half-filled and almost half-filled 
cases in the sti^ong-coupling regime. The conjecture leads to another conjecture: If |n — 1| <^ 
i} I [DU"^), the Gutzwiller band is at the center of the Hubbard gap between the upper and lower 
Hubbard bands; and if |n — 1| » t^/ (Df/^), the Gutzwiller band is at the bottom of the upper 
Hubbard band or the top of the lower Hubbard band, depending on whether n > 1 or n < 1. 

4.4. Metallic conductivity 

If ipi — +00, the density of states at the chemical potential, p{0), is vanishing. It is interesting 
to examine whether or not the conductivity is vanishing in such a case. 
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Magnetic impurities are introduced into the Hubbard model: 

n' = -^Jl[cr'^-' ■S';)dld^^„ (76) 

i(T(T' 

where S[ is an impurity spin at the ith unit cell. An ensemble is considered for the exchange- 
interaction constant J/. It is assumed that is positive, or zero, or negative, and that it is com- 
pletely random from unit cell to unit cell and from sample to sample: 



= 0, plj;)) = 6^j\J'\\ (( J;jj4» =0, . . . , (77) 

where ((•••)) stands for the ensemble average. If once the ensemble average is taken, the trans- 
lational symmetry is restored in the averaged system. Thus, the self-energy due to impurity 
scatterings is diagonal with respect to the wave number. In this subsection, it is assumed that 
\J'i\ ^ ksTK- In the Bom approximation, the self-energy is given by solving self-consistently 

'^^{iei) = 4>''^S'iS' + I)\jfjy2i-g,iiei,k), (78a) 

g^{iei,k) = - — — — — - — . .— .. . , (78b) 

lei + fi* - ^(fc) - jk{i£i) - {l/(t>i)I^a{m) 

where S' is the magnitude of the impurity spins, Efj{i£i) is the ensemble-averaged self-energy, 
which corresponds to r{iei) introduced in ^ and g^{iei,k) is the ensemble-averaged one of 
grj{iei,k), which is defined by Eg. (l60l). Since Ws = (^s/<^i = 0(1), as shown in Eg. (l29l). 
{l/4>i)Sa{i£i) is of the zeroth order in l/0i. According to Eg. (I34l ). 0i/[7r^(O)] is of the zeroth 
order in l/c^i , so that 7k {isi ) is also of the zeroth order in I /(pi. 

According to the Kubo formula,!^ the electrical conductivity is given by 

(Txx (w) = — [K^x (io + iO)- K^AO)], (79) 
where K^x (w + iO) is the retarded one of 

K^i^i) = dre'^^^ (e^^ie-^^i) = A_^n,^iioji), (80) 



where is the first or x component of the current operator defined by 

j = ^jik)hk, 3{k) = -^^E{k), nfc. = i^e*MH.-fi.)4d^.^, (81) 

fccr ij 

and Uxxii^i) is defined by 

n^^{ii^i) = -^^s\n{k^a)s\u{p^a) I dre"^''^ {e'^^hkae~"^^hp^') , (82) 

kp aa' •'^ 

where = ki and px = Pi- The current vertex jx{k) has to be consistently renormalized with 
Ua {i^i ) and AS^^^^ {iei , fe) in order to satisfy the Ward relation.ES The vertex correction due to 
impurity scatterings vanishes in the Born approximation, while the ladder vertex of the first order 
in the Js{iuji, q) has to be considered. If Egs. (1471 ) and (l53l) are used, it follows that 

1 2'Kxx{ii^l) 

ll + {?,J/2D)W^TTxx{i^i) 



nxxiioji) = xxy_u ^g3-^ 



T^xxiiuJi) = —'^sm^{kxa)g^{iei, k)g^{iei +iuji,k). (84) 

nk 
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According to Eq. dTPl ). the (xJ^-linear term of Kxx{i^i) contributes to the static conductivity axx{0), 
so that the a;;-Iinear term of Uxxii'^i) contributes to it. It follows that 



1 



dUxxjiu^i) _ 2. 

According to Eq. (IC-61 ). it follows that 

l + {3J/2D)W^7Txx{0)\ =t/t*. 



dTTxx{iU}l) 



d{iuji) 



If Eqs. (1851) and (1861) are used, the static conductivity is simply given by 

e2 (2t*)2 

Sxx{^) = lim — ['Kxx{'^ + «0) - vr^x(O)] 



2h 



E 



sm 



r+oo 

/ de 






1 — oo 


de 





Imc/^(e + iO,fc) 



(85) 



(86) 



(87) 



(88) 



If 2 siv?{kxa) = 1 — cos(2A;^a) is used and the term that includes cos(2A;^a) is ignored, then 



Sxx{0) = - dE ME) 

TT 



+ CX1 

de 

oo 



fie) 



de 



Im- 



e-E- 7K(e + iO) - {l/4>i)^Ae + iO) 



(89) 



If T ^ Tk or Tn < T ^ Tk, depending on whether Tn does not exists or it exists, and if 



2tk 



de 



\ df{e) 


[-7K(e + iO 


de 


-=/ 


' + 00 

de 

-oo 


\ df{e)^ 


2ts J. 


de 



n 



[-{l/4>i)lmL^{e + iO)], 



{kBT)\ 



(90) 



(91) 



are approximately used for — 7K(e+«0) and —{l/4>i)lmU^{e+iO), respectively, and if the energy 
dependences of p{E) and ( l/0i) Reiser (e + iO) are ignored, then it follows that 



(TxxiO) 



e2 8|t*P 



MO) 



Da^-^l/TK) + (I/ts)- 



(92) 



Since \t*\ = 0{\J\/D) and <^ip(0) = 0(l/|t*|), the conductivity is of the zeroth order in 
1/01. Thus, it does not vanish even in the limit of (pi — >• +oo, although /9(0)|t| — as (^i — +oo. 
In a clean system, H/ts =0. If no symmetry is broken even at T = K in the clean system, then 
1 /tk — >• as T — K, so that the conductivity diverges as T — K. The RVB electron liquid 
shows a metallic conductivity at least under the Bom approximation^ although it is almost a spin 
liquid in the sense that the density of states at the chemical potential, p{0), is almost vanishing. 



*' If impurity scatterings are rigorously treated, even if they are weak, the conductivity of the RVB electron liquid 
in one and two dimensions has to be vanishing as T — ^ K because of the Anderson localization.'^ 
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§5. Adiabatic Continuation 

In the Heisenberg limit of U /\t\ — )• +00, with J = —At^/U being kept constant, the half- 
filled Hubbard model is reduced to 

■Hs = \L{ed- ii)+V- 
where — ^ = —U /2, the impurity term considered in ^4.4l is included, and 

^=n(i-^w)("it+^4), ^i=^E'^'"''4^^-'- ^94) 

i acr' 

Because of the projection operator V, empty and double occupancies are completely excluded. 
Since Sj's satisfy the commutation relation for spin within the constrained Hilbert subspace where 
no empty nor double occupancy is allowed, Eq. ( |93l) is none other than the Heisenberg model. 
It is easy to confirm that 

[^5, rii^ + nii\ = 0, \H, + n^] / 0, (95) 

for any ith site. The local gauge symmetry exists in the Heisenberg model, but it does not in 
the Hubbard model; and the conductivity of the Heisenberg model is absolutely zero, but the 
conductivity of the Hubbard model with no impurity can be divergent at T = OK even in the 
Heisenberg limit, unless no symmetry is broken even at T = K, as studied in ^4.41 This situation 
for the Heisenberg and Hubbard models is similar to the situation for the s-d and Anderson model. 
In the s-d limit, the Anderson model is reduced to the s-d model; the local gauge symmetry exists 
in the s-d model, but it does not exist in the Anderson model. 

The local gauge symmetry is a peculiar symmetry. The local gauge symmetry cannot be spon- 
taneously broken.'*''' In the reduction of the Hubbard and Anderson models into the Heisenberg 
and s-d models, the local gauge symmetry is not spontaneously restored but i?, forced to be restored 
by the constraint of the Hilbert space within the subspace where no empty nor double occupancy 
is allowed. One of the most peculiar features of the local gauge symmetry is that the difference 
of it is irrelevant to the adiabatic continuation between the local electron liquid in the Anderson 
model and the local spin liquid in the s-d model. 

The strength of magnetic impurities can be used as an adiabatic parameter. Here, it is assumed 
that —oo<J'-< +00 and < | J'|^ < +00. Clean and dirty limits are defined by the limit of 

I J' I — )• and the limit of | J'| — )• +00, respectively. In the dirty-limit Hubbard model, an 
electron is localized almost within a unit cell, so that the local gauge symmetry is almost restored 
and the conductivity is almost zero. Therefore, it is certain that every physical property of the 
dirty-Umit Hubbard model in the Heisenberg limit is the same as that of the dirty-limit Heisenberg 
model. Thus, the electron state in the dirty-limit Hubbard model in the Heisenberg limit and the 
spin state in the dirty-limit Heisenberg model are adiabatically connected to each other. 

According to the scaling theory for the Anderson localization,'*' ' there is no critical point 
between metallic and insulating phases, or between itinerant and localized states, or between the 
clean and dirty limits; and there is no lower limit of the metallic conductivity nor no minimum 
metallic conductivity. Therefore, the RVB electron liquid in the clean-limit Hubbard model is 
adiabatically connected to the electron state in the dirty-limit Hubbard model. It is obvious that the 
RVB spin liquid in the clean-limit Heisenberg model is adiabatically connected to the spin state in 
the dirty-limit Heisenberg model. Thus, the RVB electron liquid in the clean-limit Hubbard model 
and the RVB spin liquid in the clean-hmit Heisenberg model are adiabatically connected to each 



1 J 



V, 



(93) 
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other. The difference of the local gauge symmetry is also irrelevant to the adiabatic continuation 
between the RVB electron liquid in the Hubbard model and the RVB spin liquid in the Heisenberg 
model. 

§6. Discussion 

The Hubbard model in one dimension is particular, because no symmetry can be brokenP^) 
and because the Bethe-ansatz exact solution was given.lS) Since the Bethe-ansatz solution is for the 
canonical ensemble, the electron number is an integer. According to the Bethe-ansatz solution, 
the ground state is not the Mott insulator at least for any non-half-filling, i.e., for any N ^ L and 
even for A'^ = L zt 1. In the grand canonical ensemble, on the other hand, the averaged electron 
number {J\f) can be an irrational. If a critical 5Nc is defined such that the ground state is a metal 
for I (TV) - L| > 5Nc while it is the Mott insulator for | (TV) - L\ < 5Nc, it is anticipated that 
the critical 5Nc, if it exists, is within the range of < 5Nc < 1. In the thermodynamic limit, 
the electron density is defined by n = lim^^+oo (TV) /L. At least if 6Nc/L — )• as L — )• +oo, 
the ground state is a metal for any n / 1, even for n = 1 it 0+, in the thermodynamic limit. 
The Hubbard model with n = 1 it 0^ is none other than the half-filled Hubbard model. The 
ground state of the half-filled Hubbard model is the Tomonaga-Luttinger electron liquid^^ and 
is also the RVB electron liquid. If 6Nc/L — )• as L — )• +oo, effects of 6Nc/L being different 
from zero can be ignored in the thermodynamic limit, so that the RVB electron liquid in the half- 
filled Hubbard model with n = 1 it 0+ is adiabatically connected to the RVB spin liquid in the 
Heisenberg model. 

It is assumed that 6Nc/L — as L — +oo, and that the chemical potential n is one such that 
n = 1 ± 0+ for it, and that U/\t\ > 1 or U/\t\ +oo with J = -Af /U being kept constant. If 
the conjecture of Eq. (1741 ) is true or not, is almost infinitesimal but still nonzero. Since l/4>i 
can be defined, /cbTk is nonzero: /cbTk oc \ J\ and ^bTk = 0(| J|). Therefore, the ground state is 
a metal, provided that no impurity exists and the Anderson localization does not occur, as studied 
in ^4.41 According to the Fermi-surface sum rule,'^''^ the Fermi surface or point is given by 

A:f = ±[^/(2a)](l + 5n), (96) 

where 5n = n — 1 = ±0+. If no multisite term but the RVB self-energy is considered, the 
single-particle excitation spectrum is given by 

^{k) = cjJcos{ka), (97) 

where cj ~ 1. The fi*, which is defined by Eq. (I38cl ). can be determined by the Fermi-surface 
sum rule: fi* = cj\J\ sin(7r(5n/2). Since 6n = ±0^, = ±0+. Then, a pair excitation spectrum 
u){q) is given by 

uj{q) = cj\J\ {cos[(/c + q)a] — cos(A;a)} , (98) 
where k and q are restricted to cos[(A; + q)a] — fi* > and cos(A;a) — /u* < 0. Then, 

cj\Jsm{qa)\ < uj{q) < 2cj\Jsm{qa/2)\. (99) 

On the other hand, the spectral weight of the Gutzwiller band is almost vanishing because l/(f)i is 
almost infinitesimal. Furthermore, if U/\t\ S> 1 or U/\t\ — +oo and if 6n = ±0+, low-energy 
quantum charge fluctuations are much suppressed and vanishing. Thus, the spectral weight in the 
charge channel of the pair excitation uj{q) is vanishingly small, i.e., the pair excitation u}{q) is 
almost a spin excitation. This spin-excitation spectrum in the Hubbard model in one dimension 
resembles that in the Heisenberg model in one dimension.!^ 
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On the basis of the criterion that if a gap opens in the spectrum of adding or removing a 
single electron in the canonical ensemble then the ground state has to be an insulator, as well as 
the exact result of the Bethe-ansatz solution that the gap opens for any nonzero U, Lieb and Wu 
proposed that the half-filled ground state is the Mott insulator for any nonzero On the other 
hand, according to the Bethe-ansatz solution, the residual entropy per unit cell is zero or vanishing 
as L — )• +CXD. According to the previous paper,^'*' if a complete gap opens in the self-consistent 
p{e), the residual entropy is In 2 per unit cell; even if the RVB mechanism is considered, this 
conclusion does not change. On the basis of these two results on the residual entropy, it can also 
be proposed that the half-filled ground state cannot be an insulator with a complete gap open. The 
criterion used by Lieb and Wu in order to conclude that the ground state is the Mott insulator is not 
a sufficient condition for the ground state being an insulator but a necessary condition for it; e.g., a 
similar gap opens in a metallic fine particle because of the long-range Coulomb interaction. Thus, 
it is desirable to critically reexamine whether the half-filled ground state of the Hubbard model in 
one dimension is the Mott insulator or a metal. If it is the Mott insulator, the critical Uc, which 
is defined on the basis of Eq. ([T6l ). is zero for D = 1 and {M) = L, or n = 1; and it has to be 
determined how large or small the critical 5Nc, which is defined above, and \miL^j^^{5Nc/ L) 
are. If it is a metal, the argument above for n = 1 it 0+ is valid for the exactly half-filled case of 
n = 1. 

On the basis of the adiabatic continuation, it is anticipated that physical properties resemble 
each other between the Hubbard and Heisenberg models even in one dimension, except for phys- 
ical properties that are related to the itineracy of electrons. It has already been proposed that the 
spin liquid in one dimension can be described as the Tomonaga-Luttinger spin liquid."*^ i, 46 1,471 
existence of charge fluctuations is not a necessary condition for the Fermi or Tomonaga-Luttinger 
liquid, i.e., a normal or anomalous Fermi liquid to be stabilized, as discussed in ^ According to 
the study of the present paper, this proposal is none other than the proposal that the Tomonaga- 
Luttinger electron and spin liquids are also the RVB electron and spin liquids, respectively, and 
that they are adiabatically connected to each other. The proposal that the ground state of the 
Heisenberg model is the Tomonaga-Luttinger spin liquid never contradicts the proposal that the 
exactly half-filled ground state of the Hubbard model is the Tomonaga-Luttinger electron liquid 
and not the Mott insulator. 

The half-filled Hubbard model on the square lattice is also particular, because no symmetry 
can be broken at a nonzero and because pQ {e) diverges logarithmically as e — )• 0. It follows 
that 

1/X.(0,q;r) = [1/Xs(0;r)] - (l/4)/,(0,q;r) 

= feBTK(T){l-0[[T/TK(T)]2]} -(1/4)7,(0, q;r), (100) 

where the Fermi-liquid relation studied in ^3.2l is used. In the temperature range defined by Eq. (|7]l, 
but with 5T being sufficiently large, the logarithmic divergence of 4>ip{e) as e — )• is suppressed 
by the imaginary part of the self-energy. Then, Tk{T) ~ |J|/(2A;b) or Tk{T) = 0[| J|/(2A;b)] . 
On the other hand, if T is very low such that T ^ | J| / (2A;b ), the increase or divercence of 4>ip{e) 
as e ^ is substantial, so that Tk (T) ^ K and 1/xs (0; T) ^ as T ^ kET 

Since l/xs(0; T) — )• as T — )• OK and the superexchange interaction Js(0, q) is maximum 
at q = Q, where Q = (±1, ibl)(7r/a), the ground state un the square lattice is the Neel state 
whose ordering wave number is Q or the RVB electron liquid very close to the Neel state. If T is 



*' Since the vanishing Tk{T) as T — 5- K is simply because of the divergence of (^ip(e), it never means that the 
condensation energy of the RVB electron liquid is also vanishing as T — >■ K. The condensation energy is 0(| J|/2) 
per unit cell. 
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sufficiently lower than | J|/(2/cb). it is plausible that such a low-T phase is in the critical region. 
If the low-T phase is really in the critical region of antiferromagnetism, it is anticipated that an 
anomaly appears in the uniform susceptibility, as discussed below. It follows that 

1/X.(0, 0; T) = l/xs(0, Q; T) + (1/4) [A, + Aq{T) + Ar{T)\ , (101a) 
Z\, = J,(0,Q)- J,(0,0), Z\Q(r) = jQ(0,Q;r)- jQ(0,0;r), (101b) 
Ar{T) = ylr(0, Q; T) - Ar{Q, 0; T). (101c) 

The Js{0,q) is antiferromagnetic, i.e., Js(0, Q) > and Js(0, 0) < 0; and As = 4| J|, which is 
almost independent of T. According to previous papers,n3'E).l22l jf the density of states has a sharp 
peak at one of the band edges and the chemical potential is in the vicinity of the peak position, 
Jq(0, q; T) is strongly ferromagnetic, i.e., Jq(0, 0; T) is positively large and Jq(0, 0; T) increases 
as T — OK; if the Fermi surface shows a sharp nesting, it is strongly antiferromagnetic, i.e., 
Jq (0) Qn ; TV where Qn is a nesting wave number, is positively large and Jq (0, Qn ; T) increases 
as T — )• Ko In the half-filled Hubbard model on the square lattice, the Fermi surface shows a 
sharp nesting for Q = (±1, ibl)(7r/a), so that Jq{0,Q;T) is positively large at a sufficiently 
low T and Jq{0,Q;T) increases as T — )• K; and the density of states has a logarithmic peak at 
the chemical potential, so that Jq(0, 0; T) is also positive at a sufficiently low T and Jq(0, 0; T) 
increases as T — OK. Since the peak of the density of states is at the band center and since 
the nesting effect is larger or stronger than the logarithmic-peak effect, Aq{T) > and the T 
dependence of Jq(0, Q; T) is much stronger than that of Jq(0, 0; T). Thus, the T dependence 
of Aq{T) is large; Aq{T) increases as T — OK. The Neel temperature Tn cannot be nonzero 
because of critical fluctuations or 71(0, q; T), which means that the T dependence of yl(0, q; T) 
is large. In general, the q dependence of the mode-mode coupling term is small. Then, the q 
dependence of yl(0, q; T) is small, and Ar{T) is therefore small, and the absolute, not relative, T 
dependence of the small A r {T) is also small. In the critical region, in general, the T dependence 
of l/xs(0, Q\ T) deviates from the Curie-Weiss T dependence. If the low-T phase is really in the 
critical region, l/xs(0, Q; T) ~ has to be satisfied in it. If the T dependence Ar{T) is really 
small, therefore, it is anticipated on the basis of Eq. (IIOII ) that the T dependence of 1/xs (0, 0; T) 
resembles the T dependence of Aq{T) in the critical region. Since Aq{T) increases as T — )• K, 
it is anticipated that l/x;s(0, 0; T) increases as T — K. It is interesting to examine whether or 
not l/xs(0, 0; T) increases or Xs(0, 0; T) decreases as T — OK in the half-filled Hubbard model 
on the square lattice. If the decrease of Xs(0,0;T) as T — )• OK is true, it corresponds to the 
suppression of the static and uniform susceptibility as T — )■ K in the Heisenberg model on the 
square lattice.'^'ES) 

It is straightforward to extend the study in this paper to the triangular lattice. If T > K, no 
symmetry can be broken.^^' The electron state atOK<T<TKorOK<T< | J|/(2A:b) in the 
Hubbard model is a frustrated electron liquid in the sense that no symmetry is broken in it, and the 
spin state at K < T <C | J|/(2A;b) in the Heisenberg model is the RVB spin liquid proposed by 
Fazekas and Anderson.'^ We propose that the frustrated electron liquid in the Hubbard model is 
none other than the RVB electron liquid, and that the RVB electron and spin liquids in the Hubbard 
and Heisenberg model on the triangular lattice are adiabatically connected to each other. 

In three dimensions and higher, it is possible that Tn ^ \J\/{Dk-Q) because of frustration 
and quasi-one or quasi-two dimensionality in the Hubbard and Heisenberg models, and the devia- 

*' Since A{Q,q;T) is of higher order in 1/D, its T dependence cannot be responsible for the Curie-Weiss T 
dependence in sufficiently high dimensions. Thus, only the T dependence of Jq(0, 0; T) or Jq(0, Qn; T) can be 
responsible for the Curie-Weiss T dependence of itinerant-electron magnetism in sufficiently high dimensions; on the 
other hand, only the T dependence of Xs (0; T) can be responsible for the Curie- Weiss T dependence of local-moment 
magnetism in sufficiently high dimensions, as discussed in ^ 
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tion from the half-filling in the Hubbard model, and so on. It is interesting to study how magnetic 
properties at Tn < T ^ \J\/{Dk^) resemble each other between an electron liquid in the Hub- 
bard model in the strong-coupling regime, which is the RVB electron liquid, and a spin liquid in 
the Heisenberg model, which is the RVB spin liquid. 

\f D ^ +00 and n = 1, /cbTk oc \J\/D and A{Q,q) 0. If T/Tr +00, 
•^(3(0> q) — ^ and yl(0, q) — ^ 0. If T » Tk and n = 1, the static susceptibility of the Anderson 
model is given by 



X.(0) = l/{k^T) - O [T^/{k^T^)\ . (102) 
Thus, if D ^ +00 and n = 1, the static susceptibility of the Hubbard model is given by 



^'^^^'^^ l-(l/4)J,(0,q)xs(0) fcBr-(l/4)J,(0,q)- ""'^^ 

This Xs(0;Q) agrees with the susceptibility in the mean-field approximation for the Heisenberg 
model. The Neel temperature is given by 



Tn = J.(0, Q)/(4A:b) = I J|/(2fcB), (103b) 

where Q = (±1, ±1, • • • , ibl)(7r/a). This Tn also agrees with the Neel temperature in the mean- 
field approximation for the Heisenberg model. Since the mean-field approximation becomes rig- 
orous in the limit of D — )• +00, these agreements are quite reasonable. 

If Tn ^ Tk then itinerant-electron magnetism appears at T < Tn, while if Tn ^ Tk 
then local-moment magnetism appears at T < Tn. Thus, itinerant-electron magnetism and local- 
moment magnetism are characterized by Tn ^ Tk and Tn ^ Tk, respectively, and they are 
adiabatically connected to each other. If Tn ^ Tk, a paramagnetic phase at T > Tn is none other 
than the paramagnetic phase of the insulator proposed by Slater,^'' ' or the antiferromagnetic type 
of the Mott insulator. The entropy of the paramagnetic phase is as large as In 2 per unit cell, and 
the static susceptibility of it obeys Eq. (11031 ). at least approximately. In particular, if T> — +00 
and n = 1, then Tk/Tn — >• 0, as discussed above. Thus, magnetism for D — +00 and n = 1 
is a typical one of local-moment magnetism; and the electron state for D — )• +00 and n = 1 is a 
typical one of the antiferromagnetic type of the Mott insulator. It is obvious that the local-moment 
states in the Hubbard and Heisenberg models are adiabatically connected to each other. 

The RVB electron liquid studied in this paper is none other than a normal state in order to 
study possible low-temperature ordered phases such as the Neel state of itinerant-electron mag- 
netism, which is of the zeroth order in 1/D, an anisotropic superconducting state, which is of 
higher order in 1/D, and so on in the strong coupling regime defined by U/\t\ » 1 and in the 
half-filled or almost half-filled case. It is plausible that the normal state proposed by Anderson!^ 
is none other than the RVB electron liquid studied in this paper. The study of this paper con- 
firms the relevance of theory of high-temperature superconductivity in cuprate oxides based on 
the Kondo-lattice theory.l^'IS) If the half-filled ground state in one dimension is the Mott insula- 
tor, as was proposed by Lieb and Wu,'^ it is anticipated that the RVB electron liquid studied in 
this paper can also be used as a normal state in order to study the Mott insulator in one dimension 
as a possible low-temperature phase. It is interesting to elucidate what effect is responsible for a 
complete gap to open in the Mott insulator in one dimension; within a preliminary study, however, 
there is no evidence that the self-energy can be so anomalous that it can make a complete gap 
open. 
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§7. Conclusion 

The resonating-valence-bond (RVB) electron liquid in the half-filled or almost half-filled 
Hubbard model in the strong-coupling regime is studied by the Kondo-lattice theory. Physical 
properties of the Hubbard model, such as the self-energy, and the polarization function, and the 
three-point vertex function, are decomposed into their single-site and multisite properties. Every 
single-site property can be mapped to its corresponding local property of the Anderson model. On 
the basis of the mapping to the Anderson model, the Kondo temperature Tk or fcBTk is defined 
as the energy scale of local quantum spin fluctuations in the Hubbard model. The superexchange 
interaction can also be derived by field theory. If the onsite repulsion is U and the transfer integral 
between nearest neighbors is —t/\/D, where D is the dimensionality, the exchange interaction 
constant between nearest neighbors is J/D, where J = —4t^/U; this J agrees with the one given 
by the conventional derivation. The Fock-type self-energy, which is of the first order in the su- 
perexchange interaction, is none other than the RVB self-energy. Because of the RVB self-energy, 
/cbTk is as large as \J\/D. If D is sufficiently small such that no critical temperature Tc exists or, 
even if Tc exists, Tc <^ Tk, then the RVB electron liquid is stabilized by the RVB self-energy at 
a sufficiently low T such that T < Tk or Tc < T -C Tk. If U/\t\ is finite or unless the filling of 
electrons is exactly half, A;bTk is nonzero and finite. Thus, the RVB electron liquid, which is char- 
acterized by ^bTk > 0, is stabler than any phase characterized by /cbTk = 0. In the Heisenberg 
limit, the density of states at the chemical potential is almost vanishing, so that the RVB electron 
liquid is almost a spin liquid or a quasi spin liquid. However, the quasi spin liquid shows a metallic 
conductivity at a sufficiently low temperature, provided that impurity scatterings are sufficiently 
weak and they are treated in the Bom approximation. 

According to the previous studies on the Kondo effect, the local electron liquid in the An- 
derson model and the local spin liquid in the s-d model are adiabatically connected to each other, 
although the local gauge symmetry does not exists in the Anderson model while it exists in the 
s-d model. According to the scaling theory for the Anderson localization, if no symmetry break- 
ing nor restoration occurs in a metal-insulator transition, there can be no critical point between 
metallic and insulating phases. This fact means that the metallic and insulating phases are adia- 
batically connected to each other, even if the conductivity of the metalUc phase is divergent and 
that of the insulating phase is zero. On the basis of these previous studies and the study in the 
present paper, it is proposed that the RVB electron liquid in the Hubbard model and the RVB spin 
liquid in the Heisenberg model are adiabatically connected to each other, even if the local gauge 
symmetry does not exist in the Hubbard model and the conductivity of the RVB electron liquid 
is metallic while the local gauge symmetry exists in the Heisenberg model and the conductivity 
of the RVB spin liquid is zero. The difference of the local gauge symmetry is irrelevant to the 
adiabatic continuation. 
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On the other hand, AZ!^{e + iO, k) — > as e — )- ±00; e.g., AS^ -^(e + iO, fc) — ;> as e — )• ±00, 
as shown in Eq. (|57] ) in ^4.11 Then, it is straightforward to show that 



Appendix A 

Sum Rule for the Self-Consistent /\(e) 



An analytic function F{e + iQ) is defined by 

F{e + iO) = [e + ^ - erf - U„{e + iO)] - l/R^{e + iO). 



(A-1) 



According to the mapping condition of Eq. (I13bl ). it immediately follows that 

ImF(e + iO) = A{e). 



(A-2) 
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Since F(e + iO) is analytic in the upper-half complex plane, according to Eqs. (IA-21) and (IA-31 ). 

/+00 
(ie4(e) = 27rt^. (A-4) 
-oo 

This is none other than the sum rule for the self-consistent A{e). 

Appendix B 

Theoretical Constraints for the Self-Consistent l/tpi 

B . 1 . Lower limit of the self-consistent 

If U/\t\ » 1 and n ~ 1, the density of states p{e) has a three-peak structure with the 
Gutzwiller band between the upper and lower Hubbard bands; their band centers are e ~ and 
e ~ ±(C//2), resepctively. In general, if p{e) = — (l/7r)Imi?o-(e + iO) has a peak at an e, 
IleRa{s + iO) becomes zero at the e or a little different e from the e. Thus, Rei?(j(e + iO) = 
for e's such as e = eo — and e = e± ~ ±C//2. If Rei?CT(e + iO) = 0, Eq. ([TS] ) becomes simple: 

7rp{e)A{e) = 1 - TTp{e)\lmU„{e + iO)\. (B-1) 

If n = 1 or n = 1 it 0"*", the band center of the Gutzwiller band is at the chemical poten- 
tial; thus, eo = or eo = ±0+. According to Eq. ([28]), |ImZ'^(+iO)| = 0{(f>ikBT^ /Tk); and 
according to Eq. ([32]), p{0) = p{0) = O [l/{(f)ikBTK)] ■ Thus, according to Eq. (lB4T) . 

7r/5(0)4(0) = 1-0 [(T/Tk)'] . (B-2) 
If n = 1 orn = 1 it 0+, Eq. (I66l ) is satisfied rather than Eq. (|65] ). According to Eqs. ( [69l ) and (IB-21 ). 

^(0) = O[(0it2)/(D;7)]. (B-3) 

Since ^bTk = O , as shown in Eq. dVOll, Eq. (IB3] ) is consistent with Eq. dMll- If (pi 

is large, A{0) is large and A{e) has a peak at e = 0. Since the peak width is 0(/i;B?k) or 
0[t^/ (i^f/)] and the peak height is given by (IB-31 ). it follows that 

/ deA{e) = 0[f/{DU)] x 0[{iit^) / {DU)] = 0[iit^ / [DUf]. (B-4) 

According to the sum rule of Eq. (IA-41 ). Eq. (IB-41 ) has to be smaller than 27rt^. Thus, 

(^1 < {l/c!^){DUf/t\ 1/01 > c;tV(^f^)', (B-5) 

where c'^ > and = 0(1). Even if n = 1 orn = 1 ± 0+, there is a theoretical upper limit for 

the self-consistent 0i or a theoretical lower limit for the inverse of it, l/(/>i . 

According to Gutzwiller's theory,'"'' there is another theoretical lower limit for the self- 
consistent 1/(^1, as shown in Eq. (1721 ). Thus, at least 

1/01 > max[(c'//(^f^)',cg|n- 1|], (B-6) 

have to be satisfied. If the RVB mechanism is considered, and if U is finite, and even if n = 1, the 
self-consistent <j)i cannot be divergent; and l/4>i cannot be zero. Thus, if n = 1 or n 1 it 0+, the 
critical Uc, which is defined in ^3.1.11 exists and it is infinite; and if n / 1, it does not exist. 
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If U /\t\ S> 1 and n ~ 1, the peak height and the bandwidth of the upper and lower Hubbard 
band are p(e±) = 0{l/\t\) and Wh = 0{\t\), respectively. Since |lmZ'^(e± + iO)| = 0{\t\), 
A{e±) = 0{\t\). Then, it follows that 

«±+2|t| 

/ deA{e) ~ A{e±)WYi = 0(1^). (B-7) 

Je±-2\t\ 

This is consistent with the sum rule of Eq. (IA-41) . 

B.2. Asymptotic behavior of the self-consistent l/</>i 

The hybridization energy of the mapped Anderson model is given by 

A{e) = lY}yk\^^[e+~p-Ec{k)], (B-8) 
^ k 

where L is the number of unit cells, Vk is the hybridization matrix between localized and conduc- 
tion electrons, and Ec{k) is the dispersion relation of a conduction electron. The Fermi surface 
is defined hy ft = Ecikp)- Thus, Eq. ([T6l ) is a sufficient condition for the existence of the Fermi 
surface. 

Since the A{e) defined in terms of Vfc and Ec{k) is crucial for the Kondo effect, as shown in 
Eq. (fTTI) . it can be assumed without the loss of generality that Vk is constant such that Vk = V . 
Then 

A{e)=^\V\^~p,{e), ~p,{e) = iY,6[e+~^-E,{k)\, (B-9) 

^ k 

where /5c (e) is the density of states of the conduction band. According to the sum rule of Eq. (IA-41) . 
it follows that n\V\'^ = 2t^. 

The half-filled Anderson model, in which n = lorn = libO"'", can be mapped to the s-d 
model in the s-d or Heisenberg limit of U /\t\ — >• +oo, with J = —At^/U being kept constant; the 
exchange interaction constant is given by 

Js-d = -4|y|Vf/ = -(8/7r)(tVf^), (B-10) 

and the density of states of the conduction band is equal to the p{e) in the Anderson model. The 
dimensionless coupling constant, which is defined by 

m = Js-dPcie) = -[4/(7r?7)]Z(e), (B-11) 

is relevant for the Kondo effect in the s-d model. 

If g{e) is constant as a function of e, for example, in the most-divergent approximation,"^ the 
Kondo temperature is given by 

^bTk = VFe-i/l3(o)l, (B-12) 

where W is the half of the conduction bandwidth. If ^(0) > 0, k^TK is nonzero. Since the energy 
dependence of g{e) is the same as that of A{e), as shown in Eq. (IB-Ill ), in the mapped s-d model, 
the energy dependence of ^(e) has to be seriously considered. According to the scaling theory 
for the s-d model,^''^' high-energy processes substantially renormalize fixed-point or eventual 
low-energy properties but they can cause no symmetry breaking; and the eventual low-energy 
properties play a crucial role in the quenching of the localized spin by the Kondo effect. Thus, 
whether the eventual /cbTk is zero or nonzero depends on whether the bare ^(0) is zero or nonzero. 
If ^(0) > 0, the eventual /cb^k is nonzero; and if ^(0) = 0, the eventual /^bTk is zero. 
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If Eqs. (IB-31 ) and (IB llb are used, it follows that 

\m\=0[4>itV{DU^)]- (B-13) 

There are three possibilities fo the asymptotic behaviors of \g{0)\ and l/i^i in the Heisenberg 
limit: 

lim \g{0)\=0, lim l/(^i > ^VPf/^), (B-14a) 

(7/|t|-++oo U/\t\^+oo 

0< lim |5(0)|<+oo, lim l/4>i = 0\t'^/(DU'^)], (B-14b) 

f//|t|-i>+oo (7/|t|^+oo 

lim \g{0)\=+oo, lim l/<^i < tVPf^^)- (B-14c) 

!7/|t|->-+oo U/\t\^+oo 

Equation (IB14al ) is inconsistent with the self-consistent k^Tx being nonzero, and Eq. (IB14cl ) is 
inconsistent with Eq. (IB-51 ). Thus, Eq. (IB14bl ) has to be satisfied in a self-consistent solution for 

n = 1 and n = 1 it 0"*". 

Appendix C 

Proof of the Equality ofEq. 



There is a useful relation between Sd, which is defined by Eq. (|58] ). and the static -Kxxi^), 

which is defined by Eq. (l84l) . as studied below. In the presence of magnetic impurities, they are 
described as 

1 X /' + 00 

Sd = ryZ^oik) def{e)ling^{e + iO,k) 

= -V Vcos(feia) / de/(e)Im5,(e + iO,fc), (C-1) 

vr.x(O) = ^ sin2(A:ia) J def{e) [lm.g,{e + iO, k)] [^eg^{e + zO, fe)] . (C-2) 



Equation (IC ll ) is also given, in the integration form, by 

= ——Vn dki--- dkDCOs{kia) def{e)lmg^{e + iO,k). (C-3) 

^(^TTj J-7r/a J-Tv/a J -oo 



By the partial integration of Eq. (IC-31 ) with respect to fci, it follows that 

-.D />+7r/a />+7r/a /■+00 



Sn = 2t* / d/ci---/ dkDsin^ikia) def{e)[lmg^{e + iO,k)][Reg^{e + iO,k)]. 

(C-4) 



-n/a J —n/a 

This is also given, in the sum form, by 



2 /■+00 
SD = 2t*—Y,^\T?{kia) de/(e)[lm5,(e + iO,fe)][Re5,(e + iO,fc)]. (C-5) 

This is simply 2t*iTxx{^), so that 
Equation (l86l) is derived from Eq. (IC-61) . 



